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ABSTRACT 



Context. The majority of the intermediate mass, pre-main-sequence Herbig Ae/Be stars reside in binary systems. As these systems 
are young, their properties may contain an imprint of the star formation process at intermediate masses (2-15 Mq). However, these 
systems are generally spatially unresolved, making it difficult to probe their circumstellar environment to search for manifestations of 
their formation process, such as accretion disks. 

Aims. Here we investigate the formation mechanism of Herbig Ae/Be (HAe/Be) binary systems by studying the relative orientation 
of their binary orbits and circumstellar disks. 

Methods. We present linear spectropolarimetric observations of HAe/Be stars over the Ho- line, which are used to determine the 
orientation of their circumstellar disks. In conjunction with data from the literature, we obtain a sample of 20 binaries with known 
disk position angles (PAs). We subsequently compare our disk PA data to a model to investigate whether HAe/Be binary systems and 
their disks are co-planar Moreover, in the light of a relatively recent suggestion that some HAe/Be star spectropolarimetric signatures 
may not necessarily be related to circumstellar disks, we re-assess the relationship between spectropolarimetric signatures and disk 
PAs. We do this by comparing spectropolarimetric and high spatial resolution observations of young stellar objects (both HAe/Be and 
T Tauri stars). 

Results. We find that spectropolarimetric observations of pre-main-sequence stars do indeed trace circumstellar disks. This finding 
is significant above the 3cr level. In addition, our data are entirely consistent with the situation in which HAe/Be binary systems and 
circumstellar disks are co-planar, while random orientations can be rejected at the 2.2<t level. 

Conclusions. The conclusive alignment (at more than 3cr) between the disk PAs derived from linear spectropolarimetry and high 
spatial resolution observations indicates that linear spectropolarimetry traces disks. This in turn allows us to conclude that the orbital 
planes of HAe/Be binary systems and the disks around the primaries are likely to be co-planar, which is consistent with the notion 
that these systems form via monolithic collapse and subsequent disk fragmentation. 



Key words. Techniques: polarimetric - 
Ae/Be 



Stars: emission-line - formation - binaries - pre-main-sequence - variables: T Tauri, Herbig 



1. Introduction 

Herbig Ae/Be stars are i ntermediate mass, 2-15 Mq, pre-main- 
uence stars (see e.g. Herbig 19601: Hillenbrand et a 



seq 



L; ■ IWaters & Waelkenslll998HHernandez et alJl2004 . Thev are the 
' most massive objects to experience an optically visible pre- 
main-sequence (PMS) phase. Therefore, they provide a key op- 
^ portunity to study the early evolution of stars more massive than 
5^ \ the lower mass T Tauri stars at these wavelengths. Herbig Ae/Be 
(HAe/Be) stars span the transition between low and high mass 
objects. As a result, it has been postulated that there is a change 
in accretion mechanism within the HAe/Be mass r ange (see e.g. 
IVink etai]|2002l: [Eisner et al.ll2004tlMonnieil2005l) . If this is the 
case, the more massive Herbig Ae/Be stars may form in the same 
fashion as the most massive stars. Consequently, the study of 
Herbig Ae/Be stars may be used to constrain the process of mas- 
sive star formation, which is still not fully understood (see e.g. 
IZinnecker & Yorkell2007h . 

iBaines et al.l (|2006') used spectroastrometry to detect unre- 
solved HAe/Be binary systems and they compared the posi- 
tion angle (PA) of the systems to the orientation of their cir- 
cumprimary disks. They reported that the binary and circumpri- 



mary disk PAs of the systems in their sample are preferentially 
aligned. This indicates that the circumstellar disks lie in the same 
plane as the binary orbit. In turn, this suggests that these systems 
formed via a scenario that features fragmentation, as opposed to 
alternatives involving the capture of binary companions. 

This is consistent with one model of massiv e star formation. 
Using 3D radiation-hydrodynamic simulations, iKmmholz et al.l 
show that stars of even 40/V/q can form via monolithic 
collapse and disk accretion. A prediction of these models is that 
disk fragmentation leads to binary systems. The resultant bi- 
nary system is aligned with the original fragmented disk struc- 
ture. Therefore, because of angular momentum considerations, 
the binary system is also aligned with the inner accretion disk 
that reaches onto the primary star's surface. In addition, the sec- 
ondary, at a distance of approximately 1000 AU, has a relatively 
high mass. It has also been proposed that massive stars form 
via competitive accretion. In this scenario, a star gains mate- 
rial from its parental cluster. In the process, binary systems are 
formed via stellar capture. These binaries start as solar mass bi- 
naries at separations of the order of 1000 AU and evolve into 
close high mass binaries at approximately 1 AU dBonnell & Batd 
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I2005h . much closer than in the monolithic accretion scenario. 
Also, as these systems are formed via stellar capture, the binary 
orbital planes a nd circumprimary disks d o not necessarily lie in 
the same plane dMoeckel & Ballvll2007h . unlike the monolithic 
collapse and disk fragmentation scenario. Therefore, the proper- 
ties of binary systems can be used to constrain the star formation 
mechanism at high masses. 



The work of Baines et al. I (l2006h was extended by 
IWheelwright et all ( 1201 ol) who used spectroastrometry to deter- 
mine the mass ratios of HAe/Be binary systems. They found that 
the mass-ratio distribution of HAe/Be binary systems favours 
comparable masses, further indicating that these sy stems formed 
via di sk fragmentation. However, the sample of iBaines et al.l 
(l2006h which initially suggested this was small (6 objects), while 
a comparison with model simulations of co-planar binary or- 
bits and disks was lacking. Therefore, the relative alignment of 
HAe/Be binary orbits and circ umstellar disks is st ill uncertain. 
We aim to extend the work of iBaines et al. I (l2006h and address 
this uncertainty by enlarging the sample of binary HAe/Be stars 
with known disk PAs. 

To achieve this, we require a large sample of HAe/Be bi- 
nary systems with known disk orientations. We generate a sam- 
ple of 20 such systems using linear spectropolarimetry and high 
spatial resolution data from the literature. Spectropolarimetry 
offers a unique opportunity to probe the circumstellar envi- 
ronment of young stellar objects on scales small enough to 
study accretion disks. If a star is surrounded by a hot disk, 
free electrons in the disk can polarise the light of the cen- 
tral star Emission-line photons emanating from the disk do not 
pass through the polarising medium, unlike the continuum light. 
Consequently, emission lines are depolarised with respect to the 
continuum. Therefore, a depolarisation signature over an emis- 
sion line can be used to i nfer the presence of a small-scale, 
otherwise undetected disk ('Clarke &McLeanl 1 1974t iPoeckertI 
[l975; Poeckert & Marlborough 1976). Alternatively, if accretion 
shocks produce emission lines close to the star, the presence 
of a disk may be manifest b y scattering pola risation of a com- 
pact source of Une emission (IVink et alJl2002L hereafter V2002). 
Furthermore, spectropolarimetric signatu res can also be used to 
constrain the geometry of such disks (e.g. lVink et al. 2005b). 

Linear spectropolarimetry is well established as a tech- 
nique to study spatiall y unresolved d isks around HAe/Be stars 
(see e.g. lOudmaiier & Drew! Il999t IVink et al.1 l2002l l2005at 
iMottram et al.ll2007l) . However, the use of s pectropolarimetry to 
study disks has recently been revisited bv iHarrington & KuhnI 
(l2007h . The authors report observations of spectropolarimetric 
signatures across absorption components of Ha emission lines. 
It is suggested such signature s are due to a proc ess other than 
scattering in a disk. Instead, iKuhn et aP (l2007h propose that 
optical pumping and absorption in winds may be responsible. 
Optical pumping, i.e. a preferential population of different mag- 
netic sub-states of the lower level of an atomic transition, may 
be caused by an anisotropic radiation field. If a gas is optically 
pumped, its opacity will depend upon the orientation of the elec- 
tric field of incident radiation. As a result, the absorption of light 
in such gas can resul t in linear polarisation effects across ab- 
sorption features ('see lKuhn et al.ll2007l l2010l) . We note that all 
spectropolarimetric signatures require a flattened, asymmetric 
geometry otherwise the polarisation vectors cancel. This is in- 
dependent of the polarising mechanism. Nonetheless, if the po- 
larisation signatures of HAe/Be stars originate in significantly 
asymmetric winds, this does not automatically demand an in- 
ward coplanar extension of any large-scale circumstellar disk. 



Here we investigate whether there is a clear relationship between 
the angle of polarisation and disk orientation which, at minimum 
signals alignment between the unresolved and resolved spatial 
scales. Such an outcome is required if indeed the bulk of the 
observed linear polarisation of these objects is due to a circum- 
stellar disk. 

This paper is structured as follows. In Section[2]we evaluate 
the alignment of HAe/Be binary systems and circumstellar disks 
using linear spectropolarimetry over Ha to determine the PAs 
of disks in HAe/Be binary systems. To justify the use of spec- 
tropolarimetry to trace circumstellar disks, we reassess the rela- 
tionship between spectropolarimetric signatures and disk PAs in 
Section[3] Finally, we discuss the results and conclude the paper 
in Section|4] 

2. The relative orientation of HAe/Be binary 
systems and circumstellar disks. 

The key to determining whether HAe/Be binary planes and cir- 
cumstellar disks are preferentially aligned is to use a large sam- 
ple of binary systems with known disk PAs. There are several 
studies of HAe/Be star binary systems (see e.g. Baines et al] 
120061: [Wheelwright et alJl2()T(l . However, there are few mea- 
surements of the orientation of the circumstellar disks in these 
systems. Linear spectropolarimetry is the favoured technique to 
probe the orientation of HAe/Be star disks since a sample of ~20 
can be obtained in a few nights of observing (see e.g. V2002). 
To maximise our sample of disk PAs, we have conducted new 
Ha spectropolarimetric observations; taken linear spectropolari- 
metric results from previous work; and assimilated additional 
constraints on HAe/Be star disk orientations from the literature. 
We first present the additional spectropolarimetric observations 
we undertook (Section lTTI ). We then use the assembled sample 
of disk PAs, in conjunction with a simple model (described in 
Section 123] ). to evaluate the hypothesis that circumstellar disks 
in HAe/Be binary systems lie in the orbital plane (Section l24l i. 

2.1. Linear spectropolarimetric observations 

We present linear spectropolarimetric observations of HAe/Be 
stars conducted in the R band and centred on Ha. Here we de- 
scribe the observing procedure and the data reduction steps be- 
fore presenting the observed spectropolarimetric signatures. 

The sample wa s predominately chosen from the catalogue 
of lThe et al.1 (11994 . Targets were selected based on their visual 
magnitude (V < 12-13). This is because spectropolarimetry re- 
quires high signal-to-noise-ratio (SNR) data and thus only bright 
targets can be observed. Most of the targets were binary Herbig 
Ae/Be stars that had not been observed with spectr opolarimetry 

[ revio usly, and these were pr imarily drawn from IBaines et al.l 
2006h and iThomas et al.l (|2007). Several objects with existing 
spectropolarimetric data were observed to check the results are 
consistent with previous observations. 

The linear spectropolarimetric data were obtained using the 
WilHam Herschel Telescope (WHT) from 08-11-2008 to 10-11- 
2008. Clouds were present for the majority of the three nights, 
preventing observation for some, but not all, of the time. The 
seeing was typically fair (~1"), although on occasions it became 
poor (2-2.5")- The observations were conducted with the ISIS 
spectrograph which was equipped with polarising optics com- 
prising of a calcite block and a rotating half-wave plate. The 
R1200R grating was used and the central wavelength was set 
to 6560 A. Several slit widths were used, ranging from 1 to 1.8 
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arcsec, and the minimum spectral resolution was found to be 
-4250, or 70 km s-'. 

The calcite block separated the incident light into two per- 
pendicular rays: the ordinary (O) and the extraordinary (E) rays. 
Each observation comprised of both the O and the E ray spec- 
trum of the science target, and a corresponding set of sky spec- 
tra. The polarisation atPAs ofO°, 22.5°, 45° and 67.5° was mea- 
sured by rotating the half-wave plate. Multiple polarisation stan- 
dard stars were observed to characterise the instrumental polari- 
sation and calibrate the polarisation angle. A log of the observa- 
tions is presented in Table [1] 

Data reduction was conducted using the Image Reduction 
and Analysis Facility (irafQ in conjunction with routines writ- 
ten in Interactive Data Language (idl). The data reduction pro- 
cess for each observation consisted of trimming, bias subtrac- 
tion, flat-field division and cosmic -ray removal. Following the 
above, the target O and E spectra, and those of the sky if they 
were present, were extracted from each frame. Sky spectra, 
which were not always detected, were typically a few percent 
of the stellar spectra. We note that polarisation signatures over 
Ha, which are the focus of the paper, are unaffected by con- 
taminant sky polarisation. Therefore, the sky polarisation has no 
influence on our final results. Wavelength calibration was per- 
formed using CuNe and CuAr arc spectra, which were obtained 
periodically during the observing run. 

Once the O and E ray spectra had been extracted, the Stokes 
parameters for each data set were calculated using a routine writ- 
ten in IDL. The method used is that outlined in the ISIS polarisa- 
tion manual by Jaap Tinbergen and Rene Rutter0. For each set 
of polarisation data, i.e. data obtained with the half-wave plate 
at 0° and 45° or 22.5° and 67.5°, the ratio of the O and E rays 
in each spectrum was calculated. To obtain the degree of polar- 
isation, the data obtained at a given PA were averaged and the 
following equations were used: 



^0,45°/'^E,45° 



q ■ 



R-\ 
R + \ 



(1) 



(2) 



where q - Qjl, I is the total flux input and /o.ana and 
^E.ang are the fluxes of the O and E rays at the stated half-wave 
plate PAs. This procedure was repeated for the other set of po- 
larisation data, i.e. data obtained with half-wave-plate PAs of 
22.5° and 67.5°, to calculate u. 

To calculate the total polarisation and the polarisation angle 
the data were combined using the following equations: 



4 



-tan 



(3) 



(4) 



where P represents the total polarisation, and 9 is the polari- 
sation angle. 

Instrumental polarisation was not corrected for The standard 
observations indicate that the instrumental polarisation is of the 



' http://iraf.noao.edu/, see lTodvl l ll993h 
^ http://www.ing.iac.es/Astronomy/observing/manuals 
/html_manuals/wht_instr/isispol/isispol.html 



order 0. 1 per cent, and is not the dominant source of continuum 
polarisation. Interstellar polarisation is not corrected for either, 
as such corre ctions are typ ically subject to significant uncertain- 
ties (see e.g. lJensen et al.ll2004 ). Contaminant polarisation sim- 
ply adds a wavelength independent vector to the Stokes Q and U 
parameters. Therefore, plotting the spectrally dispersed q against 
u allows the intrinsic angle of polarisation, and hence the polar- 
ising media, to be estabhshed. 



2.1 .1 . Spectropolarimetric signatures 

We present an example of the spectropolarimetric signatures ob- 
served in Figure [T] Data around Ho' for all the stars in the sam- 
ple are presented in AppendixlAlin Figure lATl while continuum 
polarisation measurements are included in Table [1] The signa- 
tures of the objects previously observed with spectro polarimetry 
are generally consistent with publ ished results (e.g. IVink et al.l 
'2001 12005 aL iMottram et"aDl2007h . This provides an important 
check on the data reduction process. In general, the data are 
of slightly inferior quaUty to previous observations. This is at- 
tributed to the poor weather conditions throughout the observ- 
ing. Consequently a coarser binning is used than is typical for 
such data. Five objects exhibit a change in both linear polar- 
isation and the polarisation angle over Ha (HD 179218, HK 
Ori, MWC 1080, V586 Ori & MWC 147). MWC 1080 and 
MWC 147 were also observed bv lMottram et alj ( l2007l hereafter 
M2007). While the data are broadly consistent with the previous 
observations, the line effects are not obvious in the QU diagram. 
This is probably a result of the coarse binning. Less coarse bin- 
ning does not reveal any signatures as the scatter increases con- 
siderably. Therefore, the results of M2007 are used rather than 
these new data. Of the three remaining objects that exhibit line 
effects, HK Ori and V586 Ori exhibit an excursion in QU space 
and HD 179218 exhibits a clump of data surrounding the contin- 
uum value, with a slight extension along the U axis. 

To arrive at the disk position angle from these spectropo- 
larimetric signatures, one should have some knowledge of the 
polarising mechanism. In general, if the polarimetric signature 
is the result of simple depolarisation, the polarisation vector in 
QU space should be measured from the line to the continuum. In 
case of intrinsic line polarisation, the reverse is true, effectively 
resulting in a difference of 90° from the former mechanism. To 
differentiate between the two, the width of the line-effect is often 
used as a proxy (see V2002). The intrinsic polarisation angles of 
HK Ori and V586 Ori are determined via linear fits to their QU 
excursions, assuming the signature is due to depolarisation (as 
the width of the spectropola rimetric signature s is comparable to 
the emission line width, see IVink et al.l 120021) . The polarisation 
angle of HD 179218 is calculated assuming the excursion is only 
in the U direction and that the signature is due to intrinsic polar- 
isation (since the signature is narrower than the emission line). 
The uncertainties in the polarisation angles calculated are ap- 

o 

proximately 10 . 



The continuum polarisation present ed in Table [His broadl y 
consistent with literature values (see e.g. lMaheswar et aT]l2002h . 
However, there are exceptions. For example, we note that the 
continuum polarisation of HK Ori and MWC 758 differs from 
literature values. This may be due to intrinsic variability as 
both objects are known to exhibit variab le polarisation (see 
iBaines etal ] l2004l:lBeskrovnava et alJ[T999l) . 
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Table 1. The log of observations. The continuum polarisation was measured in the wavelength region 6520-6600 A, excluding 
changes over the Ha line. The uncertainty in the continuum polarisation is typically 0. 1 per cent and the uncertainty in the continuum 
polarisation angle is of the order of 1°. The polarisation angles are in the equatorial frame. 
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RA 


Dec 


Spec. Type 


V 


T 
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Seeing 


Date 


P 
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(J2000) 
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(mags) 


(mins) 
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(%) 


n 


XY Per 


HD 275877 
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9.4 


120.0 
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50 
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Fig. 1. An example of the spectropolarimetric signatures of the sample. The top half of the figure presents the spectropolarimetric 
PA, the percentage polarisation, and the Stokes intensity spectra centred upon Ha. The data are binned to a constant polarisation 
error, which is stated in the plots. The solid red line is the un-binned line profile. The QU diagrams of the signatures are displayed 
in the lower half of the figure. The solid lines mark the direction of the intrinsic polarisation. 



2.2. Literature data 

To supplement the spectropolarimetrically observed sample, in- 
trinsic polarisation angles (determined using spectropolarime- 
try) and binary PAs were taken from the literature. The resultant 
sample of objects for which both spectropolarimetric and binary 
PAs are available is presented in Table |2] As the original studies 
from which the data are taken were not primarily concerned with 
disk orientations, we re-assessed all the polarisation angles taken 
from the literature to ensure they are calculated consistently. The 
QU data associated with the PAs listed in Table|2]were employed 
to determine approximate angles. If these differed by more than 
90° from the literature values, the reported angles were rotated 
by 90° (see the discussion above on calculating disk PAs from 
spectropolarimetry). Table |2]presents the resultant sample. 



The sample was then further supplemented by adding disk 
PAs determined from direct imaging and multi-baseline interfer- 
ometry. The additional sample is also presented in Table |2] 

2.3. Modelling co-planar binary systems and disks 

In order to compare disk and binary PAs in a meaningful way, 
a model is required to predict the difference in th e two angles 
expected for various scenarios. iBaines et al ] (120061) simplv com- 
pare disk and binary PAs. However, binary PAs do not necessar- 
ily relate to the binary orbital plane. For example, even when a 
binary system and its circumprimary disk lie in the same plane, 
if the system is seen face-on, the binary PA is unrelated to the PA 
of the disk. In the case of more edge-on systems, the binary PA 
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Table 2. Binary systems for which a measurement of the intrinsic polarisation angle (from linear spectropolarimetry) or a direct 
constraint on the orientation of the circumprimary disk and the PA of the binary system is available. Column 3 denotes the spectral 
type of the system primary, taken from SIMBAD unless otherwise stated. Column 4 lists the binary PA, column 5 lists the intrinsic 
polarisation angle and column 6 contains the disk PA. Given the typical errors, the angles are presented to the closest degree. 



Object 


Alt. Name 


Type 


Bin. PA 

n 


Pol. PA 

(°) 


Disk PA 

(°) 


Only spectropolarimetry 










MWC 166 


HD 53367 


BO 


298' 


46^ 




HD 58647 


BD -13° 3008 


B9 


115^ 


20" 




MWC 158 


HD 50138 


B9 


30^^ 


135'* 




MWC 120 


HD 37806 


A2 


34' 


90" 




V586 On 


HD 37258 


A2 


217' 


81^ 




TOri 


MWC 763 


A3 


107' 


20" 




HK Ori 


MWC 497 


A4 


47' 


169^ 




Spectropolarimetry and imaging/interferometry 






HD 200775 


MWC 361 


B2 


164« 


93^ 


7' 


MWC 147 


V700 Mon 


B6 


82' 


168^ 


80'* 


HD 45677 


FS CMa 


B2'' 


150^ 


164"' 


77" 


BD +40° 4124 


V1685 Cyg 


B3 


175'2 


36- 


110'^ 


MWC 1080 


V628 Cas 


BO 


269' 


752 


55" 


CQTau 


HD 36910 


F3 


56'2 


20" 


120'" 


HD 179218 


MWC 614 


AO 


141'2 


-45= 


23'5 


Only imaging/interferometry 










MWC 758 


HD36112 


A5 


311'2 




128'^ 


V892 Tau 


HBC 373 


B8 


23 '2 




53'« 


RMon 


MWC 151 


BO 


287" 




~80'« 


MWC 297 


NZ Ser 


BO 


313" 




I6520 


HR 5999 


V856 Sco 


A7 


-HP' 




-25^2 


HD 101412 


PDS 57 


B9.5 


226'2 




38'5 



( 20 10|), 2 : M2007, 3: iBaines et al.l (l2006h. 4: IVink etalJ ( |2005ah. 5: These data, 
(,20091), 8: iKi-aus e t all ll2008h. 9: ICidale et al.l 7200100: jPatel etal] ||2006|). 11 



References. 1: From the data presented in [Wh eelwright et al 
see Section 12.1.1! 6: [Pirzkal et alj (Il997h 7~lokamoto et al 
iMonnier et all (l2006h 12: iThomas et al.l (l2007h. 13: lEisner'et al 



VI— ^ (l2004^■ 1 4: iDoucet et al. (2006), 15: Fede le et al.l (l2008h. 16 

20081), 17: IWeigelt et all ( l2002h . 18: iFuente et al.l i2006l) . 19: IVinketal.1 ( l2005ci) . 20: iManoi et alj (i2007l) . 21: IStecklum et alT 
iPreibisch et al.l ( l2006^ . 



Monnier et al.l 
T995I1 



22: 



is likely to be aligned with the disk PA. Nonetheless, at certain 
phases of the orbit, the binary PA will be quite different to the 
circumprimary disk PA. It is only in the extreme case of an edge- 
on, co-planar system that the disk and binary PAs are constantly 
aligned. Here we employ a new model to predict the average 
distribution in the difference between disk and binary PA when 
circumstellar disks and binary orbits lie in the same plane. 

The model is characterised by a random orbital phase; in- 
clination; semi-major axis and PA of the line about which the 
system is inclined. The eccentricity of the system was kept con- 
stant. Although a constant eccentricity is not truly representative 
of the eccentricity di stribution of PMS binary systems (see e.g. 
iGoodwin et al.ll2007h . neither is a random eccentricity distribu- 
tion. Since the eccentricity can affect the results, the input ec- 
centricity is left as a constant and is treated as a free parameter. 
This allows us to fit the observed distribution in differences be- 
tween binary and disk PAs. The masses of the components were 
kept constant at 6 & 1 Mq. The component masses have little 
influence on the final distribution and are mentioned only for 
completeness. 

We used 10000 random systems to determine the distribu- 
tion in the difference between the disk and binary PAs (APA). 
Spectropolarimetry is insensitive to face on systems as the pro- 
jected polarisation vectors cancel one-another out. Therefore, 
systems with low inclinations are discarded. The value of the 
cut-off inclination does affect the final distribution. However, so 
long as large values (e.g.> 50°) are not chosen, the differences 
between the final distributions are relatively small. Since we also 



use imaging observations which are less sensitive to inclination, 
a low cut-off value (10°) is used. The imaging and spectropolari- 
metric samples have different cut-off values. However, neither 
sample is large enough to be used in isolation, and thus a single 
cut-off value is required. Changing the cut-off value has a similar 
effect to changing the eccentricity. Since we use the eccentricity 
as a free parameter, the use of a single cut-off value for the sam- 
ple of both imaging and spectropolarimetric data will not prevent 
a fit to the data. By varying the eccentricity to fit the data, the un- 
certainty in the cut-off value is essentially incorporated into the 
best-fitting eccentricity as a systematic uncertainty. 

Figure |2] presents a typical distribution. We are only con- 
cerned with the magnitude of the difference between disk and 
binary PAs. Therefore, we convert all differences between these 
angles to be in the range - 90°. For example, a difference of 
170° is 10° away from alignment, and is thus classified as a dif- 
ference of 10°. An alignment between disk and binary PAs cor- 
responds to an offset of 0°. It can be seen that the distribution 
tends towards intrinsic alignment, and appears noticeably differ- 
ent to the distribution expected if the two angles are not related. 
Specifically, the distribution does not exhibit a direct correlation 
between binary and disk PA, but does demonstrate an excess of 
aligned angles over the random distribution. This is to be ex- 
pected based on the previous discussion. To reiterate, co-planar 
disks and orbits are more likely to be observed to be aligned than 
not, since at high inclinations the disk and binary PAs will be the 
same. This is not the case for non-coplanar systems. Therefore, 
the co-planar distribution exhibits a preference for an alignment 
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40 60 
APA (") 

Fig. 2. The expected distribution in the difference between po- 
sition angles of disks and binary systems in a co-planar model 
(blue solid line) with an eccentricity of e = 0.5. The black dotted 
line represents the expected distribution due to a random orien- 
tation of binary and disk PAs. 



between disk and binary position angles compared to the random 
distribution. 



We note that comparing this model to observations assumes 
that the orientation of the circumstellar disks in the sample is 
primordial. This need not be the case. Even if the circumstellar 
disks in HAe/Be binary systems originally lie in the plane of the 
binary orbits, gravitational interactions with another star have 
the potential to alter their orientation (see e.g. the discussion in 
iBate et al.ll200(]|) . This would weaken an intrinsic correlation be- 
tween binary and disk PA. Therefore, any observed correlation 
between disk and binary PA may be more significant than it ini- 
tially appears. However, quantifying this is beyond the scope of 
the paper. 

2.4. The relationship between binary and disk position angle 

The final sample of HAe/Be binary systems with disk PAs con- 
tains 2 objects, and is thu s more than three times greater than 
that of iBaines et al ] (I2006h . The differences between disk and 
binary PAs in the sample are compared to the co-planar model 
and a random distribution. In Figure [3] we show the cumulative 
distribution of the differences in PA between the binary systems 
and those derived for the disks. Also in the figure are the distri- 
butions predicted by the co-planar model and a completely ran- 
dom association of disk and binary PAs. A first glance indicates 
that the aligned scenario provides a better fit to the data than 
the random distribution. This is corroborated by a Kolmogorov- 
Smirnov (KS) test; while the co-planar model is found to be con- 
sistent with the data (it can only be ruled out at the 0.2cr level), 
the random hypothesis is barely consistent with the data, and can 
be rejected at a level of 2.2cr. 

This 2<T level of rejection may seem low, but we note that 
both models are more similar in their cumulative distributions 
than one might naively expect. As mentioned earlier, this is be- 
cause inclined binary systems do not necessarily have the same 
observed orientation as their circumprimary disks. This may ex- 
plain why our significance, with an improved data-set, is of the 
same order as iBaines etaL I (l2006l) found with a smaller data-set 
and a more simplistic model. 
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Fig. 3. The difference in disk and binary PAs presented in Table 
|2](red dashed). This is compared to the co-planar model with a 
constant eccentricity of 0. 1 (blue solid) and a distribution gen- 
erated by assuming the two angles are randomly aligned (black 
dotted). Note that APA refers to the difference between binary 
and disk position angle (i.e. accounting for the 90° offset be- 
tween spectropolarimetric signatures and disks). 

Improving the current statistics requires further enlarging 
the sample. We find that if the data followed the model distri- 
bution exactly, differentiating between the two scenarios at a 
3cr level would require a sample of approximately 50 objects. 
Nonetheless, the current data are consistent with the co-planar 
hypothesis and furthermore, this scenario is favoured over ran- 
dom alignments. 



3. Does spectropolarimetry really trace 
circumstellar disks? 

Our use of spectropolarimetry assumes that the spectropolari- 
metric signatures of HAe/Be stars can be used to trace the ori- 
entations of their circumstellar disks. Linear polarisation has 
been used for som e time to probe the orien tation of circumstel- 
lar disk s (see e.g. Bastien & Menardlll990l). However, the pro- 
posal of [Harrington & KuhnI (l2007h and lKuhn et al.l (l2007h that 
the spectropolarimetric signatures of HAe/Be stars are due to op- 
tical pumping and not scattering in a disk implies that these sig- 
natures do not have to trace disks. We address this issue here and 
do this by comparing spectropolarimetric observations of young 
stellar objects with independent measurements of the orientation 
of their disks. To increase the sample we include observations of 
both HAe/Be and T Tauri stars. 

Table |3] presents the sample of young stellar objects, both 
low mass T Tauri objects and Herbig Ae/Be stars, for which the 
combination of spectropolarimetric observations and indepen- 
dent observations of their circumstellar disks is available. We 
first concentrate on the HAe/Be sample and then assess whether 
there is a general trend within the combined sample of young 
stellar objects. The m ajority of HAe/Be disks are thought to be 
optically thin (see e.g. iNatta et al.ll200ll) . If the observed polar- 
isation signatures are due to single scattering in these disks, the 
polarisation angles will be perpendicular to the disks' PA on the 
sky. Therefore, if the spectropolarimetric signatures of HAe/Be 
stars are due to disks, the difference between the disk and po- 
larisation angles for this sample will be 90°. In testing this hy- 
pothesis, we allow for an uncertainty in the difference between 
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polarisation and disk PAs of ~15°, which is typical for the val- 
ues used. In some cases, there is a systematic uncertainty due 
to different PAs being rep orted by different authors. F or exam- 
ple, in the case of AB Aur. [Mannings & Sargent! (Il997h report a 
disk PA of 79 while I Corder et al.l (l2005l) report various possi- 

o o 

ble angles ranging from 26 to 85 . However, such cases are not 
common, and are unlikely to influence conclusions regarding the 
whole sample. 



3.1. Herbig Ae/Be Stars 

In the case of the HAe/Be stars in Table|3] the majority of the ob- 
jects (7 out of 9) have polarisation angles approximately perpen- 
dicular to their imaged disks. This might be expected if the polar- 
isation signatures are due to single scattering in disks. Therefore, 
this appears to validate the use of spectropolarimetry to trace 
disks. Here we quantify this. We use the HAe/Be star data in 
Table[3]to test the hypothesis that the intrinsic polarisation angle 
is always perpendicular to the disk PA (within the errors). The 
cumulative distribution of the difference in disk PA and polarisa- 
tion PA is shown in Figure|4] We find that the hypothesis that the 
disk and polarisation angles are unrelated to each other and thus 
randomly oriented can be discounted at a significant level (at 
3.1cr according to the KS test). In contrast, the hypothesis that 
the spectropolarimetric signatures of the HAe/Be stars are ori- 
ented perpendicularly to their disks cannot be rejected beyond 
the Icr level, and is thus consistent with the data. Therefore, we 
find that spectropolarimetric signatures of HAe/Be stars do trace 
the orientation of their circumsteflar disks. 

We note that, although the majority of HAe/Be objects ex- 
hibit a difference in disk and polarisation angle that is close to 
90°, two objects have disk and polarisation angles that are es- 
sentially aligned. This is contrary to expectations based on sin- 
gle scat tering occurring in an optically thin disk. In their smaller 
sample, IVink et al] ( 12005 al) also note several objects where this 
is the case. They suggest that while the spectropolarimetric sig- 
natures of all young stellar objects are due to circumstellar disks, 
the angle of the resulting polarisation vector is dependent upon 
the properties of the inner disk. If the inner disk is optically thin, 
single scattering dominates and the resulting polarisation vector 
is perpendicular to the disk PA. In contrast, if the inner disk is 
optically thick, the polarisation vector is parallel to the disk PA 
due to multiple scattering. 

Many of the T Tauri stars in Table [3] also have disk and po- 
larimetric PAs that are aligned. Therefore, if this argument is 
correct, it would appear that the majority of T Tauri star disks 
are optically thick in the inner regions. 

3.2. Removing fhe dependence on optical deptli 

Although the previous test provides a statistically significant re- 
sult, it could be argued that this is dependent upon prior knowl- 
edge of the disks' optical depth. Assuming HAe/Be star disks 
are optically thin is in part justified. The vast majority of the 
HAe/Be stars considered have an offset between disk and polari- 
sation angle that is close to 90°, as expected for single scattering 
in optically thin disks. However, it could be argued that this is a 
circular argument. Moreover, in the case of the T Tauri objects, 
the disk and polarisation angles are generally aligned. A more 
appropriate hypothesis for the total sample may be a combina- 
tion of the two scenarios, i.e. that polarisation angles are related 
to the PAs of disks, but that the polarisation angles can be either 



aligned or perpendicular to the disk. In general, we do not have 
prior knowledge of the optical depth of the disks. Therefore, we 
test the relationship between polarisation and disk angle without 
making a priori assumptions. 

In the scenarios mentioned above, the difference between the 
disk and spectropolarimetric angle is either or 90°, depending 
upon the optical depth of the inner disk. Consequently, it can be 
expected that the offset from 45° to the difference between disk 
and spectropolarimetric angle (henceforth A*!*) is always 45°. 
We note that this is the case regardless of whether the signature 
is interpreted as being due to line polarisation or depolarisation, 
or whether the disk is optically thin or thick. Therefore, this test 
is even more robust than the previous which assumes optically 
thin scattering and is subject to the polarisation signature being 
interpreted correctly. Here we compare this hypothesis to the 
sample presented in Table [3] The disk and spectropolarimetric 
angles in Table [3] are used to calculate A*?. This is then com- 
pared to the hypothesis that A^P is 45° by calculating the average 
of 10000 equally sized samples in which A*!* is 45° but with an 
additional random error contribution with a maximum value of 
15°(see FigureEli. 

The hypothesis that A*? is random can be discounted at a 
significant level (above 3cr). This leaves the hypothesis that A*!* 
is always 45°, which cannot be rejected at greater than a Icr 
level and is thus consistent with the data. This confirms the ear- 
lier finding that spectropolarimetric signatures trace disks and 
demonstrates that this is not dependent on spectral type and clas- 
sification (i.e. T Tauri star or HAe/Be star). 

To summarise, we show that the spectropolarimetric data 
presented in Table[3]do appear to trace circumstellar disks. This 
would be expect ed if the polarisation i s due to scattering in these 
disks. However, [Harrington & Kulinl (l2009h claim many of the 
HAe/Be stars in our sample have polarisation signatures which 
require optical pumping and absorption in outflows (e.g. AB Aur 
and MWC 480). We show here that regardless of the polarising 
mechanism, spectropolarimetry can be employed to trace the ori- 
entation of circumstellar disks. 

4. Discussion and Conclusion 

4.1. On tlie use of spectropolarimetry to probe circumstellar 
disks 

The results presented in this paper indicate a direct correla- 
tion between the spectropolarimetric signatures of pre-main- 
sequence stars and the orientations of their circumstellar disks. 
This is significant above the 3cr level and appears indepen- 
dent of the classification of the young stellar objects. Therefore, 
we conclude that spectropolarimetric signatures of young stel- 
lar objects do indeed trace the orientation of their circumstel- 
lar disks. This is expected in the case of polarisation of stellar 
and accretion sh ock photons by disks (iMcLean & Clarkjll979l : 
IVink et aT]l2002l) . In the case of polarisation via optical pump- 
ing and absorption, the relationship between spectropolarimetric 
signatures and disks is less clear 

We note that iKuhn et aH (1201 Ol) suggest that the polarisa- 
tion signature of the Herbig Be star HD 200775 is due to op- 
tical pumping and that the signature does trace the orienta- 
tion of an imaged disk. However, the spectropolarimetric sig- 
nature of this object is observed across a double-peaked emis- 
sion line profile, and might therefore be the result of depolari- 
sation after all. Nevertheless, many Herbig Ae/Be stars, several 
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Table 3. Young stellar objects (column 1) for which spectropolarimetric observations and a direct measurement of their disk PA are 
available. The disk and adopted polarisation angles are presented in columns 4 and 5 and the difference between them is listed in 
column 6. Finally, column 7 indicates objects where the difference is close to 90° (±) or 0° (||). 



Object 


Alt. Name 


Type 


Disk PA 


Pol. PA 


APA 








(°) 


D 


(°) 


HAe/Be 












HD 200775 


MWC 361 


B2 


7' 


932 


86 _L 


MWC 147 


V700 Mon 


B6 


80' 


168- 


88 _L 


HD 45677 


PS CMa 


B2'' 


iT 


164* 


87 _L 


BD +40° 4124 


VI 685 Cyg 


B3 


110^ 


36^ 


74 _L 


MWC 1080 


V628 Cas 


BO 


55' 


752 


20 II 


CQTau 


HD 36910 


F3 


120** 


20' 


80 _L 


MWC 480 


HD 31648 


A3 


150'' 


55' 


85 _L 


AB Aur 


HD 31293 


AO 


7910 


160' 


81 _L 


HD 179218 


MWC 614 


AOlVe 


23" 


~45'2 


22 II 


T Tauri 












RYTau 


HD 283571 


F8 


62" 


163' 


79 _L 


SU Aur 


HD 282624 


G2 


127''' 


130' 


3 II 


FUOri 


HBC 186 


G3 


4715 


45' 


2 II 


GWOri 


HD 244138 


G5 


56'« 


(60)' 


4 II 


DRTau 


HBC 74 


K5 


128''' 


120' 


8 II 



References. 1 : lOkamoto et alj (l2009h . 2: M20 07. 3 : iKraus et all d2008h. 4: ICidale et al J ll2001h, 5: iMonnier et all (l2006h, 6: iPatel et all ( I2006I) . 7: 

Eisner etal. (2004), 8: Doucet et al. (2006), 9: Vink et al. (2005a), 10: Mannings & Sargent ( 1997h. 1 1 : [Fedele et al.' ( 2008|), 12: these data, 13: 
Akeson et al.. (.2003,) , 14: ,Akeson et al.,t2002.) . 15: ,Malbet et al.. (.2005,) , 16: .Mathieu et al..(.1995h . 17: lKitamura et al..(,2002ir 
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Fig. 4. The distribution in the difference between spectropolarimetrically predicted disk PA and observed disk PA for the sample 
presented in Table[3](blue dashed). This is compared to a random distribution (black short dotted). On the left we show a distribution 
where polarisation signatures are always orientated perpendicularly to circumstellar disks and on the right we present the distribution 
for a scenario where the spectropolarimetric signatures can be either perpendicular or parallel to disks (see the text for more detail). 
Both model distributions have a maximum error of 15°. In both cases, a random orientation of disk and polarisation position angles 
can be discarded at the 3cr level. 



of which are in our sample, exhibit polarisation signatures as- 
sociated with P Cygni line profiles, i.e. outflowing gas (see e.g. 
[Harrington & Kuhn 2009). Here we show that the signatures still 
appear to trace the orientation of circumstellar disks. This im- 
plies that, if the spectropolarimetric signatures are due to optical 
pumping and absorption in a wind, the wind geometry essen- 
tially mirrors that of the disk, at least in the regions where the 
polarisation occurs. 

This is partly substantiated by recent observations of 
the Hq' emission of the Herbig Ae star AB A ur with the 
CHARA array by iRousselet-Pe rraut et all (l2010l) . This object 
has been proposed to exhibit polarisation due to optical pump- 
ing since it displays polarisation a cross the P Cygni absorp - 
tion component of its Ho' emission dHarrineton & Kuhnll2007h . 



iRousselet-Perraut et a resolved the Ha emitting region 

around AB Aur and found that it could be modelled as the base 
of a wind represented by a flattened torus encompassing a cir- 
cumstellar disk. Provided the inclination and the angle between 
the wind surface and disk mid-p lane is low (20° and 35° in th e 
case of the disk-wind model of iRousselet-Perraut et a i] l2oToh . 
such a flattened torus might well appear to have a similar mor- 
phology to the disk. 

4.2. On {he alignment between binary systems and their 
circumprimary disks 

To investigate the relative alignment of HAe/Be binary systems 
and circumstellar disks we have used spectropolarimetry and 
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high spatial resolution data to determine the orientation of cir- 
cumstellar disks around the primary components of such sys- 
tems. We then combined these disk angles with binary param- 
eters to assess whether HAe/Be circumstellar disks and binary 
systems are co-planar Studies of lower mass T Tauri stars have 
found that the circumstellar disks in T Tauri star binary sys- 
tems tend to be aligned, suggesting that such systems may form 
via fragment ation (see e.g. Wolf et al. 2001; Jensen et al. 2004f 
lMoninetalJ i2006). Here we investigate whether this is also the 
case for HAe/Be systems. We note that lMaheswar et al ] (I2002h 
also compared HAe/Be binary and polarisation angles (although 
these were calculated via broadband polarimetry and thus sub- 
ject to uncertainties in the correction for the interstellar polari- 
sation). These authors find that their data are inconsistent with 
a random association of disk and binary position angles with a 
significance of 84 per cent. This is not very conclusive and the 
authors note that they do not account for projection effects and 
thus the actual correlation may be stronger. We do account for 
projection effects and demonstrate that, in principle, co-planar 
and randomly orientated disks and binaries can be differentiated. 

We show that the data are best fit with a model in which 
the binary orbit and circumprimary disk are co-planar. This is 
consistent with the suggestion that these systems formed via the 
monolithic collapse of a core and subsequent disk fragmenta- 
tion, whichjsjiowjmssive binary systems are thought to form 
(see iKrumholz et al]|2009h . However, as of yet, a random asso- 
ciation of disk and binary planes can only be excluded at a 2cr 
level. A sample of approximately 50 objects is required to reject 
the random hypothesis and thus to differentiate between the two 
scenarios at 3cr or higher. 

4.3. Final remarks 

To summarise, we find that spectropolarimetric signatures of 
young stellar objects do trace the orientation of their circum- 
stellar disks. This is independent of a specific mechanism for 
the linear polarisation. In itself, this finding is insensitive to the 
polarising mechanism as all mechanisms require some form of 
asymmetric geometry. We note that scattering in a disk appears 
a plausible polarisation mechanism as it naturally explains the 
relationship between disk and polarisation angle. Furthermore, 
assuming the polarisation of T Tauri stars is due to multiple scat- 
tering in optically thick disks, it is consistent with the observa- 
tion that T Tauri star polarisation is generally parallel to imaged 
disks while the reverse is true for most HAe/Be stars. It is not 
clear how optical pumping and polarisation via absorption can 
reproduce the different polarimetric behaviour of HAe/Be and 
T Tauri stars. Further modelling is required to investigate this 
issue. 

We conclude that our results are entirely consistent with 
the disks and orbits of HAe/Be binaries being co-planar, and 
thus with the scenario of binary formation via disk fragmen- 
tation. Further spectropolarimetric observations, e.g. provided 
by SALT, in conjunction with additional high resolution data, 
e.g. provided by NIR interferometry, are required to increase the 
sample and conclusively differentiate between aligned and ran- 
dom orientations of disks and binaries. 

References 

Akeson, R. L., Ciai'di, D., & van Belle, G. T. 2003, in SPIE Conf. Ser., ed. 

W. A. Traub, Vol. 4838, 1037-1042 
Akeson, R. L., Ciardi, D. R., van Belle, G. T., & Creech-Eakman, M. J. 2002, 

ApJ, 566, 1124 



Baines, D., Oudmaijer, R. D., Mora, A., et al. 2004, MNRAS, 353, 697 
Baines, D., Oudmaijer, R. D., Porter, J. M., & Pozzo, M. 2006, MNRAS, 367, 
737 

Bastien, P & Menard, F. 1990, ApJ, 364, 232 

Bate, M. R., Bonnell, I. A., Clai'ke, C. J., et al. 2000, MNRAS, 317, 773 
Beskrovnaya, N. G., Pogodin, M. A., Miroshnichenko, A. S., et al. 1999, A&A, 
343, 163 

Bonnell, I. A. & Bate, M. R. 2005, MNRAS, 362, 915 

Cidale, L., Zorec, J., & Tringaniello, L. 2001, A&A, 368, 160 

Clarke, D. & McLean, I. S. 1974, MNRAS, 167, 27P 

Corder, S., Eisner, J., & Sai'gent, A. 2005, ApJL, 622, L133 

Doucet, C, Pantin, E., Lagage, P O., & Dullemond, C. P 2006, A&A, 460, 117 

Eisner, J. A., Lane, B. R, Hillenbrand, L. A., Akeson, R. L., & Sargent, A. I. 

2004, ApJ, 613, 1049 
Fedele, D., van den Ancker, M. E., Acke, B., et al. 2008, A&A, 491, 809 
Fuente, A., Alonso-Albi, T., Bachiller, R., et al. 2006, ApJ, 649, L119 
Goodwin, S. P., Kroupa, P., Goodman, A., & Burkert, A. 2007, in Protostars and 

Planets V, ed. B. Reipurth, D. Jewitt, & K. Keil, 133-147 
Harrington, D. M. & Kuhn, J. R. 2007, ApJL, 667, L89 
Harrington, D. M. & Kuhn, J. R. 2009, ApJS, 180, 138 
Herbig, G. H. 1960, ApJ Supplement, 4, 337 

Hernandez, J., Calvet, N., Briceno, C, Hartmann, L., & Berlind, P. 2004, AJ, 
127, 1682 

Hillenbrand, L. A., Strom, S. E., Vrba, R J., & Keene, J. 1992, Astrophysical 
Journal, 397, 613 

Jensen, E. L. N., Mathieu, R. D., Donar, A. X., & DuUighan, A. 2004, ApJ, 600, 
789 

ICitamura, Y., Momose, M., Yokogawa, S., et al. 2002, ApJ, 581, 357 

Kraus, S., Preibisch, T., & Ohnaka, K. 2008, ApJ, 676, 490 

Krumholz, M. R., Klein, R. I., McKee, C. F, Oflfner, S. S. R., & Cunningham, 

A. J. 2009, Sci, 323, 754 
Kuhn, J. R., Berdyugina, S. V., Fluri, D. M., Harrington, D. M., & Stenflo, J. O. 

2007, ApJL, 668, L63 
Kuhn, J. R., Geiss, B., & Hanington, D. M. 2010, ArXiv e-prints 
Maheswar, G., Manoj, P, & Bhatt, H. C. 2002, A&A, 387, 1003 
Malbet, R, Lachaume, R., Berger, J., et al. 2005, A&A, 437, 627 
Mannings, V. & Sargent, A. I. 1997, ApJ, 490, 792 
Manoj, P, Ho, P T. P, Ohashi, N., et al. 2007, ApJL, 667, L187 
Mathieu, R. D., Adams, R C, Fuller, G. A., et al. 1995, AJ, 109, 2655 
McLean, I. S. & Clarke, D. 1979, MNRAS, 186, 245 
Moeckel, N. & Bally, J. 2007, ApJ, 656, 275 
Monin, J.-L., Menard, R, & Peretto, N. 2006, A&A, 446, 201 
Monnier, J. D., Berger, J.-P, Millan-Gabet, R., et al. 2006, ApJ, 647, 444 
Monnier, J. D., Tannirkulam, A., Tuthill, P G., et al. 2008, ApJL, 681, L97 
Monnier, J. D. e. a. 2005, ApJ, 624, 832 

Mottram, J. C, Vink, J. S., Oudmaijer, R. D., & Patel, M. 2007, MNRAS, 377, 
1363 

Natta, A., Prusti, T., Neri, R., et al. 2001, A&A, 371, 186 

Okamoto, Y. K., Kataza, H., Honda, M., et al. 2009, ApJ, 706, 665 

Oudmaijer, R. D. & Drew, J. E. 1999, MNRAS, 305, 166 

Patel, M., Oudmaijer, R. D., Vink, J. S., Mottram, J. C, & Davies, B. 2006, 

MNRAS, 373, 1641 
Pirzkal, N., Spillar, E. J., & Dyck, H. M. 1997, ApJ, 481, 392 
Poeckert, R. 1975, ApJ, 196, 777 

Poeckert, R. & Marlborough, J. M. 1976, ApJ, 206, 182 

Preibisch, T., Kraus, S., Driebe, T., van Boekel, R., & Weigelt, G. 2006, A&A, 
458, 235 

Rousselet-Perraut, K., Benisty, M., Mourard, D., et al. 2010, A&A, 516, LI 
Stecklum, B., Eckart, A., Henning, T., & Loewe, M. 1995, A&A, 296, 463 
The, P S., de Winter, D., & Perez, M. R. 1994, A&AS, 104, 315 
Thomas, S. J., van der Bliek, N. S., Rodgers, B., Doppmann, G., & Bouvier, J. 

2007, in lAU Symp, Vol. 240, lAU Symp., ed. W. I. Hartkopf, E. F Guinan, 

& P Harmanec, 250-253 
Tody, D. 1993, in A.S.P Conf.Ser, Vol. 52, Astronomical Data Analysis 

Softwaie and Systems II, ed. R. J. Hanisch, R. J. V. Brissenden, & J. Barnes, 

173-183 

Vink, J. S., Drew, J. E., Harries, T. J., & Oudmaijer, R. D. 2002, MNRAS, 337, 
356 

Vink, J. S., Drew, J. E., Harries, T. J., Oudmaijer, R. D., & Unruh, Y. 2005a, 

MNRAS, 359, 1049 
Vink, J. S., HaiTies, T. J., & Drew, J. E. 2005b, A&A, 430, 213 
Vink, J. S., O'Nein, P M., Els, S. G., & Drew, J. E. 2005c, A&A, 438, L21 
Waters, L. B. R M. & Waelkens, C. 1998, ARA&A, 36, 233 
Weigelt, G., Balega, Y. Y, Hofmann, K.-H., & Preibisch, T. 2002, A&A, 392, 

937 

Wheelwright, H. E., Oudmaijer, R. D., & Goodwin, S. P 2010, MNRAS, 401, 
1199 

Wolf, S., Stecklum, B., & Henning, T. 2001, in lAU Symp., Vol. 200, The 



9 



H.E. Wheelwright et al.: The relative alignment of Herbig Ae/Be binaries and disks 



Formation of Binary Stars, ed. H. Zinnecker & R. Mathieu, 295-304 
Zinnecker, H. & Yorke, H. 2007, ARA&A, 45, 481 



Appendix A: Observed spectropolarimetric 
signatures 

Here we present the spectropolarimetric signatures of the entire 
sample over Ha. Figure IaTI presents the signatures in terms of 
the total flux, the amount of polarisation and the polarisation an- 
gle over the Ho- Hne. The QU diagrams of MWC 1080 and MWC 
147, objects which display a spectropolarimetric signature that 
was not presented in the main body of the paper, are presented 
in Figure IAT2I 
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Fig. A.l. The spectropolarimetric signatures of the sample. For each object the spectropolarimetric PA, the percentage polarisation, 
and the Stokes intensity spectra are presented centred upon Ha. The data are binned to a constant polarisation error, which is stated 
in the plots. The solid red line is the un-binned line profile. 
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Fig. A.2. The QU diagrams of MWC 1080 and MWC 147. The data are binned to a constant polarisation error, which is stated in 
the plots. 
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